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Abstract. A perfect Euler cuboid is a rectangular parallelepiped with integer edges 
and integer face diagonals whose space diagonal is also integer. The problem of find- 
ing such parallelepipeds or proving their non-existence is an old unsolved mathemat- 
ical problem. The Diophantine equations of a perfect Euler cuboid have an explicit 
S3 symmetry. In this paper the cuboid equations are factorized with respect to their 
S3 symmetry in terms of multisymmetric polynomials. Some factor equations are 
calculated explicitly. 



1. Introduction. 

The search for perfect cuboids has the long history since 1719 (see [1-39]). This 
history is presented as an adventure story in [37]. Let x\, x 2 , x 3 be the edges of a 
cuboid and let di, d 2 , d 3 be its face diagonals. Then we have the equations 

{xif + (x 2 ) 2 - (d 3 ) 2 = 0, (d 3 ) 2 + (x 3 ) 2 - L 2 = 0, 

(x 2 f + (x 3 ) 2 - (dO 2 = 0, {d\f + (xi) 2 — L 2 = 0, (1.1) 

(x 3 ) 2 + ( Xl ) 2 - (d 2 ) 2 = 0, (d 2 ) 2 + (x 2 ) 2 -L 2 = 0, 

where L is the space diagonal of the cuboid. In the case of a perfect Euler cuboid the 
equations (1.1) constitute a system of Diophantine equations with respect to seven 
variables x\, x 2 , x 3 , d\, d 2 , d 3 , and L. In [40] the equations (1.1) were reduced to a 
single Diophantine equation with respect to four especially introduced parameters 
a, b, c, and u. On the base of this equation in [41] three cuboid conjectures were 
formulated. These conjectures are studied (but not yet proved) in [42-44]. 

In the present paper we apply a quite different approach to the equations (1.1). 
The equations (1.1) possess a natural S3 symmetry Indeed, the symmetric group 
S3 is composed by transformations of the set of three numbers {1,2,3}: 




(1.2) 

The transformation (1.2) is applied to the equations (1.1) as follows: 

a(xi) = x<ri, <x(di) = d^ a(L) = L. (1.3) 
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Looking at (1.1) and (1.3), one can easily see that the system of equations (1.1) in 
whole is invariant with respect to the transformations a G S3. 

The main goal of this paper is to factorize the equations (1.1) with respect to the 
S3 symmetry (1.3). We reach this goal by deriving some new equations from (1.1). 
These new equations are written in terms of the values of so-called multisymmetric 
polynomials (they generalize well-known symmetric polynomials) . 

2. Multisymmetric polynomials. 

Multisymmetric polynomials, which arc also known as vector symmetric polyno- 
mials, diagonally symmetric polynomials, McMahon polynomials etc, were initially 
studied in [45-51] (see also later publications [52-65]). Let's consider a set of vari- 
ables arranged into some m x n matrix as follows: 



M 



In 



(2.1) 



The symmetric group S n acts upon the matrix (2.1) by permuting its columns: 

a(xij) = x io j. (2.2) 

Definition 2.1. A polynomial p G Q[xn, . . . ,x mn ] is called multisymmetric if it 
is invariant with respect to the action (2.2) of the symmetric group S n , i. e. if 



p(xi 



,) =p(xu, ... ,x mn ) for all a G S n . 



Let q{x) = q(xn, ... ,x mn ) be an arbitrary polynomial of the variables com- 
posing the matrix (2.1). Then we can produce a multisymmetric polynomial by 
applying the symmetrization operator S to the polynomial q(xn, ... , x Tl 



S(q(x n , . . . ,x mn )) = — - y~] q{x 1A , . . . ,x man ). (2.3) 
n! * — ' 

(T£S„ 

Regular symmetric polynomials (see [66]) correspond to the special case m = 1 in 
the definition 2.1. Like in the case m — 1, in general case m > 1 there are elemen- 
tary symmetric polynomials. However, in this general case elementary symmetric 
polynomials are enumerated not by a single index, but by a multiindex: 

a = [ai, • • • , a m ], where a% ^ and \a\ = ot\ + . . . + a m ^ n. (2-4) 

Let's denote through x a the following monomial: 

X Xn ' ■ • ■ ' Xi ai • X2ai-\-l ' • • ■ " ^2ai+a2 ■ • ■ ^mn— a m +l ' ' ' *^mn • (^-5) 

V v ' V v ' V v ' 

OL\ OL2 (Xrn 

The variables in the product (2.5) are taken from n consecutive columns of the 
matrix (2.1). The initial group of ai of them is taken from the first row of this 
matrix, the next group of «2 of these variables is taken from the second row and so 
on. The last group of a rn variables is taken from the last m-th row of the matrix 
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(2.1). If on = 0, then the corresponding i-th group in (2.5) is empty and hence the 
variables of i-th row do not enter the monomial (2.5) at all. 

Definition 2.2. An elementary multisymmetric polynomial e a (in, . . . , x mn ) cor- 
responding to the multiindcx (2.4) is produced from the monomial (2.5) by means 
of the symmetrization operator (2.3) according to the formula 



e a( x ) = —7 S(x a ), where a). = ot,\\ ■ . . . ■ a m \. (2.6) 
a! 



Note that the ratio n\/a\ in (2.6) is always an integer number and e a (x) is the 
sum of exactly n\/a\ monomials produced from the monomial (2.5) by means of the 
permutations of variables (2.2). In the case of the trivial multiindcx = [0, ... ,0] 
the formulas (2.5) and (2.6) reduce to the following ones: 

x° = 1, e = 1. 

Like in the case of regular symmetric polynomials, there is the following funda- 
mental theorem for multisymmetric polynomials. 

Theorem 2.1. The elementary multisymmetric polynomials (2.6) with multiindi- 
ces < \a\ ^ n generate the ring of all multisymmetric polynomials, i. e. each 
multisymmetric polynomial p <E Q[xn, . . . , x mn ] can be expressed as a polynomial 
with rational coefficients through these elementary multisymmetric polynomials. 

The proof of the fundamental theorem 2.1 can be found in [51]. Unfortunately 
the elementary multisymmetric polynomials (2.6) are not algebraically independent 
over Q for m > 1 (see [60]). For this reason the expression of p as a polyno- 
mial with rational coefficients through the elementary multisymmetric polynomials 
e a {x), which is claimed by the fundamental theorem 2.1, is not unique. 

2. Multisymmetric polynomials associated with a cuboid. 

Note that the formulas (1.3) can be treated as a special case of the formulas 
(2.2). Indeed, let's compose the 2x3 matrix 



M 



Xi 


X2 


X3 


di 


di 





(3.1) 



Due to (1.3) the transformations a 6 S3 act as permutations of columns upon the 
matrix (3.1). Applying the definition 2.1 to the matrix (3.1), we get the concept 
of a multisymmetric polynomial of six variables x\, X2, x 3 and di, d 2l <^3- Now we 
calculate the elementary multisymmetric polynomials corresponding to the matrix 
(3.1). The first three of these polynomials are 

e[i, ] = xi + x 2 + x 3 , 

e [2 , ] = X1X2 +£20:3 +X3X1, (3.2) 
e [3 ,o] = xi x 2 x 3 . 



It is easy to see that the polynomials (3.2) coincide with the regular symmetric 
polynomials of the three variables x\, x 2l x 3 . The next three elementary multisym- 
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metric polynomials are similar to (3.2). They are 

e [0 ,i] = di + d 2 + d 3 , 

e [0 ,2] = d\ d 2 + d 2 d 3 + d 3 di, (3.3) 
e [0 . 3 ] = di d 2 d 3 . 

The polynomials (3.3) coincide with the regular symmetric polynomials of the three 
variables d\, d 2 , d 3 . The rest of the elementary multisymmetric polynomials are 
actually multisymmetric. They include variables from both rows of the matrix M: 

e [2 ,i] = xi x 2 d 3 + x 2 x 3 di + x 3 X\ d 2 , 

e llA] = x\ d 2 + d\ x 2 + x 2 d 3 + d 2 x 3 + x 3 d\ + d 3 xi, (3.4) 
e [1]2] = Xi d 2 d 3 + x 2 d 3 di + x 3 di d 2 . 

The polynomials (3.2), (3.3), and (3.4) constitute the complete set of elementary 
multisymmetric polynomials associated with the matrix (3.1). 

4. The first four factor equations. 

Note that the variables x\, x 2 , x 3 and d±, d 2 , d 3 in the matrix (3.1) are not in- 
dependent. They are related to each other by means of the polynomial equations 
(1.1). For this reason the elementary multisymmetric polynomials (3.2), (3.3), and 
(3.4) produced from the variables xi, x 2 , x 3 and di, d 2 , d 3 gain more algebraic rela- 
tions in addition to those present in the case of independent variables Xi, x 2 , x 3 and 
di, d 2 , d 3 (see comments to fundamental theorem 2.1). These algebraic relations 
are written as polynomial equations with coefficients in Q: 

P( e [l,0] ; e [ 2 ,0] , G[3,0] J e [0,l] ) e [0,2] , e [0,3] J e [2,l] , 6(1,1], e [l,2] , L) — 0. (4-1) 

The polynomial equations of the form (4.1) derived from (1.1) as well as those 
fulfilled identically due to (3.2), (3.3), and (3.4) are called factor equations of the 
cuboid equations (1.1) with respect to their S 3 symmetry. Our present goal is to 
reveal some of these factor equations explicitly. 

The equations (1.1) are quadratic with respect to their variables. For this reason 
it is quite likely that there are no linear relationships between multisymmetric 
polynomials (3.2), (3.3), and (3.4). As for higher order relationships, they do 
actually exist. In order to reveal them we need to consider squares, cubes, fourth 
powers etc, and various mutual products of the multisymmetric polynomials (3.2), 
(3.3), and (3.4). For the square (e [1|0] ) 2 we have 

(e[i, ]) 2 = x\ + x\ + x\ + 2 (xi x 2 + x 3 Xi + x 3 x 2 ). (4.2) 

On the other hand, from the cuboid equations (1.1) we derive 

x\+xl + xl = L 2 . (4.3) 

Applying (4.3) to (4.2) and comparing the result with (3.2), we derive 



( e [i,o]) 2 — 2e[ 2 ,o] — 0. 



(4.4) 
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The equation (4.4) is the first and the most simple factor equation produced from 
the cuboid equations (1.1). 

The polynomial e [0>1] in (3.3) is very similar to e [1]0] . For its square we have 

(e [li0] ) 2 = d\ + dl + d$ + 2 (d 1 d 2 + d 3 d 1 + d 3 d 2 ). (4.5) 

On the other hand, from the cuboid equations (1.1) we derive 

df + dl + d 3 = 2L 2 . (4.6) 

Applying (4.6) to (4.5) and comparing the result with (3.3), we derive 

(e [o , 1] ) 2 -2e [o , 21 -2L 2 = 0. (4.7) 

The equation (4.7) is the second factor equation produced from the cuboid equations 
(1.1). It is equally simple as the equation (4.4). 

In order to derive the third factor equation from the cuboid equations (1.1) we 
consider the cube (e [lj0] ) 3 and apply the formula (3.2): 



(en.o]) 3 = 4 + x\ + x 3 3 + 3*i {x\ + x%) + 

+ 3 X2 {x\ + x\) + 3 x 3 (x\ + x\) + 6 xi x 2 x 3 . 

Using the cuboid equations (1.1), we derive the following formulas: 



(4.8) 



4 


= L 2 


-dl 


4 


= L 2 Xl 


— d 2 X\ 




= L 2 


-d 2 2 , 


;r 2 


= L 2 x 2 


- dl x 2 


4 


= L 2 


-4, 


3 


= L 2 x 3 


- d 2 x 3 



(4.9) 



Substituting (4.9) into the equality (4.8), we obtain the formula 

(e [li0] ) 3 = -(xi df +x 2 d 2 + x 3 dl) + 7 L 2 (xi + x 2 + x 3 ) - 
— 3 x\ {d 2 + d s ) — 3 x 2 (d 2 + d 2 ) — 3 x 3 (df + d 2 ) + 6 x\ x 2 x 3 . 



(4.10) 



The right hand side of the formula (4.10) is a multisymmctric polynomial. For this 
reason we can apply the theorem 2.1 to it. As a result we get 



(4.11) 



( e [i,oj) 3 — 2 B[i t 2] + 6 ep i(l ] + 4 e[ 0i2 ] cp.oj — 
— 2 e[o,i] e [i,i] — e [i,o] e fo,i] ^~ ^ e [i.°i L 2 ■ 

Note that the equation (4.7) can be resolved with respect to e [0i2] : 

e l0 . 2] = ^f-L 2 . (4.12) 

Applying (4.12) to (4.11), we can write (4.11) as follows: 

2 e [i,2] + 6 C[3,o] — 2 e[ 0|1 j e[ ltl ] + e[i i0 j e 2 a j + 3 j L 2 — e 3 j = 0. (4.13) 
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The equation (4.13) is the third factor equation derived from the cuboid equations 
(1.1). It is more complicated than (4.4) and (4.7). 

There is another way for transforming the cube (e[ 10 j) 3 given by the formula 
(4.8). Indeed, we can resolve the left column of the equations (1.1) with respect to 
(xi) 2 , (x 2 ) 2 , and (X3) 2 . As a result we get 

2 d 2 2 +d 2 -d? 2 d 2 + d?-d 2 2 d? + d 2 - d 2 , A 1 „, 

x 2 = _2 1 1 j ^2 = _3 i 2 ) x 2 = _l 1 3 (4 w) 

The formulas (4.14) can be used instead of the formulas in the left column of (4.9). 
Applying these formulas to (4.8), we can get an expression analogous to (4.10) 
and then we can continue transforming it in a way similar to (4.11) and (4.12), 
expecting to get some new equation similar to (4.13). But actually we get the 
equation coinciding with (4.13). 

Now let's consider the cube (e [01] ) 3 . It is given by the following formula: 

( e[0il] ) 3 = di + dl + di + s dt (di + dt) + 

+ 3 d 2 (dl + df) + 3 d 3 [dl + dl) + 6 di d 2 da- 
Note that the equations of the left column of (1.1) can be resolved with respect to 
(di) 2 , (d 2 ) 2 , and (d^) 2 . They yield the equalities 

d 2 = x\ + x 2 , df = x\ di + x\ di, 

dl=x\+x\, d\ = x\ di + x\ di, (4.16) 

d\ = x\ + x\, d\ = x\ d?, + x\ d 3 - 

Substituting (4.16) into the equality (4.15), we obtain the formula 



(e [(M] ) 3 = 6 di x\ + 6d 2 x 2 , + 6 d 3 x\ + 6 d\ d 2 d 3 + 

+ 4di (x\ +xj) + 4d 2 (xj + x\) +4d 3 (xj+xj). 



(4.17) 



The formula (4.17) is analogous to the formula (4.10). Its right hand side is a 
multisymmetric polynomial. For this reason we can apply the theorem 2.1 and get 

( e [o,i]) 3 = 2 ep,i] + 6 £[o,3] — 2 o] e^ ij — 10 ep t o] e [o,i] 6 e [o,i] e fi,o] ■ (4-18) 
Note that the equation (4.4) can be resolved with respect to e [2 ,o]'- 

^m = \cIm~\ l2 - ( 4 - 19 ) 
Applying (4.19) to (4.18), we can write the equality (4.18) as follows: 

2 e ( 2> i] + 6 e [0 3 ] — 2 e [lj0 ] + e [o,i] e fi,o] ^ e [o,i] ^ 2 — e fo.i] = 0- (4.20) 

The equation (4.20) is the fourth factor equation derived from the cuboid equations 
(1.1). It is similar to the equation (4.13). 
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dl 


= L 2 


- T 2 


dl 


= L 2 


— x 2t 


4 


= L 2 


- T 2 





= di 


L 2 


- di x\, 




= d 2 


L 2 


- d 2 x 2 , 




= d 3 


L 2 


- d 3 x\. 



The equations of the second column in (1.1) can also be resolved with respect to 
(di) 2 , (d 2 ) 2 , and (c?3) 2 . Using them, we can write the formulas 



(4.21) 



The formulas (4.21) can be used instead of the equations (4.16). As a result we 
get another sequence of equations. However, the ultimate result appears to be 
coinciding with the equation (4.20). 

5. More factor equations. 

In the next step we consider the square (e [2 , ]) 2 . Using the formulas (3.2), we 
get the following explicit expression for this square: 

(e[ 2 ,o]) 2 = x 2 x 2 + x\ x 2 + x 2 x 2 +2 x\ x 2 x 3 + 2 x\ x 3 x\ + 2 x 2 xi x 2 . (5.1) 

In order to transform (5.1) we use the formulas (4.9). This yields 

(e [2 0] ) 2 = d\ dl + dl dl + dl d\ - 2 L 2 (d\ + d\ + d 2 3 ) - 2 (n x 2 dl + 



■ x 2 x 3 d 2 + x 3 x\ dl) + 2L 2 (x\ x 2 + x 2 x 3 + x 3 x\) + 3 L 4 



(5.2) 



The right hand side of the formula (5.2) is a multisymmctric polynomial. For this 
reason we can apply the theorem 2.1 to it. As a result we get 



2 _ 2 4 2 . 

( e [2,0] ) — — 6[0,1] 6[0,3] + g e [i,o] e [i,2] — g e [o,i] e [2,i] — 77 e 



,2 

[1.1] 



2 8 2 

[0,1 

2 



+ 2 e [o,i] e [i,i] e [i,o] + 77 e [2,o] e [ 0| 2] — 77 e [o,i] e [2,o] + 2 e [2i o] (5-3) 



— 77 sfi.o] e [o,2] + e 2 2] + 4 e [0l2 ] L 2 — 2 e 2 01] L 2 + 3 L 4 . 

Note that we can use the equation (4.20) in order to express e [0 3] through the other 
elementary multisymmctric polynomials in (4.20): 

6(0,3] = — 77 e [2,i] + 77 e (i,o] 6(1,1] + 77 tfo.i] ~ g e [o,i] e fi,o] — 77 e [0 ,i] L 2 . (5.4) 

Apart from (5.4), we apply the formulas (4.12) and (4.19) to (5.3). Then we get 

8 e[i, ] e [li2 ] — 8 e [0)l] e [2A] — 8 e [M] + 4 e 2 a] e 2 0] — 

- 4,i] - 3 4,o] + 10 4,0] ^ 2 + 4 4,i] ^ 2 + Li = °- 

The equation (5.5) is the fifth factor equation derived from the cuboid equations 
(1.1). It is more complicated than all of the previous factor equations. 

Now let's consider the other square (e [0 2 ]) 2 . Using the formulas (3.3), we get the 
following explicit expression for this square: 

(e [0 , 2] ) 2 = dl dl + dl dl + dl df + 2 df d 2 d 3 + 2 dl d 3 di + 2 dl di d 2 . (5.6) 
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In order to transform (5.6) we use the formulas (4.21). This yields 

(e [0 .2]) 2 = x 2 1 x 2 2 +x 2 2 x 2 3 + x 2 3 x 2 1 - 2 L 2 (x 2 1 +x 2 +x 2 )~2(d 1 d 2 x 2 + 

(57) 

+ d 2 d 3 x\ + d 3 di x%) + 2L 2 (d x d 2 + d 2 d 3 + d 3 d{) + 3 L 4 . 

The right hand side of the formula (5.7) is a multisymmctric polynomial. For this 
reason we can apply the theorem 2.1 to it. As a result we get 

( ,2_ 4 2 2 2 

(, e [0,2] ) — e [l,0] e [3,0] — 77 e [l,0] e [l,2] + 77 e [0,l] e [2,l] — 77 e [l,l] + 
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3 



+ 77 e [l.o] e [o,l] e [l,l] + 77 e [2,o] e [o,2] — 77 e [o,l] e [2,o] + 4 e [2.o] L + (5-8) 



2 2 9 9 2 2 4 

~ 2 e [i, ] e [o,2] + e [2 0] + 2 e [0 2] i — 2e [10] L +3X. 

Note that we can use the equation (4.13) in order to express e [3 , 0] through the other 
elementary multisymmctric polynomials in (4.13): 

1 1 1 3 1 2 1 T i 

6[ 3 ,o] — — g e[i,2] + g e [°- 1 ] e i 1 . 1 ) + g e [i,o] — g e [i,o] e [o,i] — 2 e[1 ' 01 ' ' ^ 

We apply the formulas (5.9), (4.12), and (4.19) to (5.8). As a result we get 

—8 e [lt0 j e [1)2] + 8 e [0>1] e [2il] — 8 e 2 x l] + 4 e [01] e 2 0] — 

- 4,o] - 3 4a] + 2 4c,!] ^ - 2 4,o] L 2 — 5 L 4 = 0. 

The equation (5.10) is the sixth factor equation derived from the cuboid equations 
(1.1). It is similar to the equation (5.5). 

The terms e [li2] , e [2 ,i], and e [ i jl ] are mixed in (5.5) and (5.10). We can separate 
e [12] and e [2A ] from e^ ^ in the following two equations: 



8e [li0] e [li2] — 8e [0-1] e [2il] + e [01] — ef 10] — 

- 8 e 2 0A] L 2 + 6 ef 1>0] L 2 + 3 L 4 = 0. 

4e 2 -2e 2 e 2 + e 4 + e 4 - 

^ °[1,1] [0,1] [1,0] T [0,1] ^ [1,0] 

- 6e 2 , 1] L 2 -2e 2 , o] i 2 + L 4 = 0. 



(5.11) 
(5.12) 



The above equations (5.11) and (5.12) are just linear combinations of the fifth and 
the sixth factor equations (5.5) and (5.10). 

In order to derive the next factor equation we consider the product e [2 ,o] 6[3,oj- 
Using the formulas (3.2), we get the following explicit expression for this product: 

6[2,0] e [3,o] = X\ X 2 X 3 + X 2 X 3 X x + X 3 X x X 2 . (5.13) 

In order to transform (5.13) we use the formulas (4.9). This yields 

e [2 ,o] e [3 , ] = xi d\ d\ + x 2 d 3 d\ + x 3 d[ d% - L 2 xi (d 2 + d 2 ) - 
- L 2 x 2 (d 2 + d 2 ) - L 2 x 3 (d 2 + dl) + L 4 {xx + x 2 + x 3 ). 
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The right hand side of the formula (5.14) is a multisymmetric polynomial. For this 
reason it can be expressed through elementary multisymmetric polynomials: 

e [2,0] e [3,0] = e [0,3] + 6 [1.2] C[0,2] + e [l,2] L + 

T 6[0,2] C[1,0] Ll ~~ e [l,l] e [0,l] "T e [l,0] -k • 

Now we transform (5.15) with the use of the formulas (5.9), (5.4), (4.19), and (4.12): 



_ 4 &[i t o] e fi,i] + 4 e [i,i] e [2,i] — 2 S[i t i] e f .i] ^ e [! 



[0,1] 



+ 2 e [1]2 ] e (1 0] + e [10] e [Q1] — e [10] — 2 e [12] L + (5.16) 
+ 5 e[i )0] ef 01] L 2 + 4 e [10] L 2 - 3 e [1-0] L 4 = 0. 

Then we apply (5.12) to (5.16) in order to eliminate the term with the square e 2 ^ 

4 e [ltl] e [2il] — 2 e [Ii:L] e [01] + 6 e [12] e [0>1] + 
+ 2 e [li2] e 2 _ 0] — ej^ Qj e 2 a] + e [li0] e [01] — 2 e [li2] £ 2 — (5-17) 
— e [i,o] e fo,i] + 2 ej 5 ^,-,] L 2 — 2 e [I>0 ] £ 4 = 0. 

The equation (5.17) is the seventh factor equation derived from the cuboid equations 
(1.1). Its order is higher than the order of the equations (5.5) and (5.10). 

The eighth factor equation is derived similarly. In order to derive it we consider 
the product e [0j2 ] e [0 3] . Applying the formulas (3.3) to this product, we get 

e [0 ,2] e [0 , 3 ] = d\ d 2 d> 3 + d 2 d 2 d\ + d 3 d\ d| ■ (5.18) 
In order to transform (5.18) we use the formulas (4.21). This yields 

ep.2] e [0 , 3 ] = di x\ x\ + d 2 x 3 x\ + d 3 x\ x\ — L 2 d x (x 2 + x\) — 



L 2 d 2 (x 2 3 + x\) - L 2 d 3 {x\ + x 2 2 ) + L 4 (di +d 2 + d 3 ). 



(5.19) 



The right hand side of the formula (5.19) is a multisymmetric polynomial. For this 
reason it can be expressed through elementary multisymmetric polynomials: 



'[0,2] e [0,3] — — e [l,l] e [3,0] + e [2,l] C[2,0] + e [2,l] L + 
T e [2,o] e [0,l] Ll — 6[i j0 ] i/ T 6[0,1] • 



Transforming (5.20) with the use of the formulas (5.9), (5.4), (4.19), (4.12), we get 

—4 e [0)1] e 2 ;1] + 4 e [li:L] e [1|2] — 2 e^^j e [10] + 6 e [2]l] e 2 0] + 

+ 2 ep,!] ef 0jl] + ef 0(I] e 2 10] - ef 01] + 2 e [3il] L 2 - (5.21) 
- 2 e [1;1] e^.o, L 2 + 4 e [0il] e 2 [10] L 2 + 7 ef 01] L 2 - 4 e [01] L 4 = 0. 

Then we apply (5.12) to (5.21) in order to eliminate the term with the square e 2 

4 e [ltl] e [li2] — 2 e [1A] e [1>0] + 6 e [2jl] e 2 0] + 
+ 2 e [2il] e 2 01] — e^ 01] e 2 0] + e [0 ,i] 0] + 2 e [2>1] i 2 — (5.22) 
— 2 e^j] e[ lj0 ] Z/ 2 + 2 e[ 0il j e 2 j £ 2 + ^ L 2 — 3 £[o,i] L 4 = 0. 
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The equation (5.22) is the eighth factor equation derived from the cuboid equations 
(1.1). It is similar to the equation (5.17). 

6. Concluding remarks. 

One can continue deriving factor equation more and more. In order to reasonably 
terminate this process we need some theoretical considerations. Note that the left 
hand sides of the factor equations are polynomials from the ring 

Q[ e [i,oj; e [2,oj, e^o], e [0 ,i], e [0t2 ], e [0 ,3], ep,i], e [i,2], L], (6-1) 

where e [1:0] , e [2i0] , e [3i0] , e [0 ,i], e [0 ,aj, e [0 , 3 ], e [2 ,ij, e [1A] , e [1>2 ], and L are treated as 
independent variables. If we continue deriving factor equation endlessly, their left 
hand sides would generate a certain ideal J in the ring (6.1). By means of the 
formulas (3.2), (3.3), and (3.4) this ideal J is mapped onto some certain ideal J sym 
of the ring of multisymmetric polynomials SymQ[xi, X2, X2, d\, d,2, c?3, L]. The ideal 
l sym is produced as the intersection 

I S ym = InSymQ[xi,X2,X2,di,d2,d 3 ,L], (6.2) 

where / is the ideal of the polynomial ring Q[xi, X2, X2, di, g?2, g?3, L] generated by 
the left hand sides of the cuboid equations (1.1). Calculating the intersection (6.2) 
is an algorithmically solvable computational problem. It will be considered in a 
separate paper. 
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